Quantum Friction in Nanomechanical Oscillators at Millikelvin Temperatures 
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We report low-temperature measurements of dissipation in megahertz-range, suspended, single- 
crystal nanomechanical oscillators. At millikelvin temperatures, both dissipation (inverse quality 
factor) and shift in the resonance frequency display reproducible features, similar to those observed 
in sound attenuation experiments in disordered glasses and consistent with measurements in larger 
micromechanical oscillators fabricated from single-crystal silicon. Dissipation in our single-crystal 
nanomechanical structures is dominated by internal quantum friction due to an estimated number of 
roughly 50 two-level systems, which represent both dangling bonds on the surface and bulk defects. 
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Understanding dissipation in nanomechanical oscilla- 
tors is essential to their study as well as to their ap- 
plications that are of fundamental and technical inter- 
ests 0, 0, ■ Applications of classical micro- and nano- 
mechanical systems — for instance, as single spin detec- 
tors Q and spin-torque oscillators memory elements 
6j , or as ultrasensitive sensors of biological interactions 
3 or fundamental forces such as gravity 0, Q — will 
be limited by energy dissipation, usually quantified in 
terms of the quality factor (Q). The realization of a me- 
chanical quantum system and its coherent manipulation 
will depend on the limitations imposed by decoherence 
and dissipation. There is a further motivation in ex- 
ploring mechanical and elastic effects in nanomechani- 
cal structures of mesoscopic system size; specifically, in 
multi-scale modelling of mechanical systems consisting of 
around 100 million atoms, which can now be fabricated 
and experimentally characterized comprehensively |icj . 

Central to the study of dissipation is the identifica- 
tion of the dominating dissipation mechanism and its 
scaling with system size. Some of the identifiable mech- 
anisms are energy loss due to clamping, thermoelastic 
processes, gas friction, mode anharmonicity, metal film 
friction, and coupling to dislocations and localized inter- 
nal defects. Since nanomechanical systems possess large 
surface-to-volume ratios, these mechanisms are treated 
essentially as surface effects. Of particular interest is the 
study of energy dissipation or internal friction by quan- 
tum tunneling of two-level systems that become effective 
typically below a temperature of 1 Kelvin. Despite the- 
oretical interest, experimental studies of internal friction 
in nanomechanical systems at millikelvin temperatures 
are yet to be reported. 

Klciman ct al. fill ] have measured dissipation in 
macroscopic centimeter-scale single-crystal silicon res- 
onators at millikelvin temperatures. The temperature 
dependence is similar to that in vitreos silica, which can 
be essentially understood in the framework of the so- 
called glass model of two-level systems [HE1E1 Their 
results were subsequently reinterpreted by Phillips in 



light of asymmetry distribution produced by local strains 
in crystals Greywall et al. [3] have identified an 

anomalous contribution to the low temperature dissipa- 
tion in single-crystal silicon resonators in the form of 
distinct periodic peaks. Carr et al., 01 Olkhovets et 
al., Evoy et al., and Carr and Craighead [2(| 
have characterized intrinsic dissipation in nano-elcctro- 
mechanical resonators (NEMS) in the 1-10 MHz range, 
identified with "surface and near-surface phenomena" . 
Yang et al. [ilj have studied the surface contributions to 
dissipation in ultrathin single-crystal silicon cantilevers. 
However, none of the nano-scale dissipation experiments 
were done at millikelvin temperatures. 

In this Letter, we report detailed measurements of 
frequency shift (5///) and dissipation (Q" 1 ) of a set 
of suspended megahertz-range nanomechanical oscilla- 
tors at temperatures down to 60 mK. We extract fre- 
quency shift and dissipation from the Lorentzian re- 
sponse of the resonators by driving them on resonance 
in the linear regime. The real and imaginary parts of 
the susceptibility function x( w o) = x'( w o) + *x"(^o) — 
pv 2 (2Sf // + iQ^ 1 ) are proportional to Sf/f and Q _1 
at frequencies close to the resonance frequency • The 
susceptibility function quantifies the system response due 
to its coupling to an environment, which results in a loss 
of both energy (dissipation) and quantum coherence (de- 
coherence). Further motivation for doing measurements 
on resonance is to enhance the sensitivity of system re- 
sponse to small dissipative forces. 

Both the relative shift in resonance frequency Sf/f 
and dissipation Q _1 , quantified by the full width at half 
maximum (FWHM) of the resonance peak, demonstrate 
reproducible temperature dependence. We find that the 
temperature dependence of frequency shift follows the 
expected logarithmic behaviour down to the lowest mea- 
sured temperature of 60 mK without saturation. This is 
consistent with the Phillips model of disordered crystals 
and the glass model of tunneling two-state systems. On 
the other hand, the temperature dependence of dissipa- 
tion does not show the expected T or T 3 dependence |22j . 
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FIG. 1: (Color online) a) SEM micrograph of the 7 /im 
nanobeam. The diagram shows the measurement circuit. To 
drive the beam, we pass an alternating current I(u) across the 
beam of length L in the presence of magnetic field B, to set up 
a force Fd r (u>) = I{uj)LB that is orthogonal to both. The cen- 
ter beam displacement x for the fundamental mode shape in 
the magnetic field B induces a voltage V em f(ui) = ^LBujqx(uj) 
in the gold electrode. We detect the response voltage using 
a RF network analyzer, b) The fundamental resonance mode 
of the beam at 14.586 MHz with Q = 33000. The induced 
voltage Vemf is a Lorentzian peak on top of a white noise 
background, which is due primarily to the input noise of our 
preamplifier, and c) the corresponding phase. 



Our single-crystal structures act essentially as surface 
systems due to their large surface-to- volume ratio. This 
results in a broad distribution of tunnel-splitting ener- 
gies or glass-like behavior. Furthermore, by quantitative 
comparison to the glass model, we find that the dissipa- 
tion in our nanomechanical structures is due to the inter- 
nal friction of roughly 50 two-state atoms and dangling 
bonds. However, at temperatures below 100 mK, the 
structures demonstrate excess temperature-independent 
dissipation, which cannot be explained with either the 
glass or the Phillips model fl5| . 

We fabricate our nanomechanical beams from single 
crystal Si/Si02/Si wafers using e-beam lithography and 
nanomachining. We have measured two different beams 
with lengths 6 and 7 /im (width 300 nm and thickness 
200 nm for both beams). The top surface of the beams is 
covered with a thermally evaporated 10 nm Cr underlayer 
and 80 nm Au film. We measure the structures magne- 
tomotively at the center of a 16 Tesla superconducting 
magnet in a dilution refrigerator. 

We observe the fundamental resonance modes of the 
two beams at frequencies 12.028 and 14.586 MHz. A 
typical linear response of the beam at 14.586 MHz with 
quality factor Q=33000 is shown in Fig. lb. Dissipation 
is defined as the inverse of the quality factor Q = u>/Au>. 
The displacement of the beam on resonance varies lin- 
early with the driving force according to Hooke's Law 
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FIG. 2: (Color online) a) The resonance peak, which is 
Lorentzian in the linear regime, assumes the asymmetric non- 
linear shape when the driving power is increased above -95 
dBm (trace D on the plot), b) Linear dependence of the in- 
duced power Pi n on driving power Pd r at 4 Tesla. The induced 
power corresponds to the resonance displacement x, shown 
versus the driving force on a log-log scale. The measured ef- 
fective spring constant is k e /f — 182 N/m. c) The response 
of the oscillator increases linearly with magnetic field, as ex- 
pected when the driving force F = ILB is kept constant at 
20 pN. We keep the driving force constant in order to inves- 
tigate the effect of the magnetic field on energy dissipation, 
d) Magnetic field dependence of dissipation. The data set is 
taken at 60 mK. It exhibits quadratic field dependence for 
high magnetic fields in both samples. Dissipation saturates 
in the low field region. 



x(ojo) — QF^r (u>n ) I L/ / [23| . The effective spring con- 
stant k e ff = 182 N/m is inversely proportional to the 
linear fit in Fig. 2b. 

The response V em f on resonance ui ~ loq is given by 
V em /(w Q ) = gL m f Q Idr(vo), according to the magneto- 
motive scheme. We measure the induced voltage of the 
response at constant force F^ r = ILB = 20 pN (Fig. 2c) 
to verify the linear B dependence of the response. Fig. 2d 
shows magnetic field dependence of dissipation in two 
samples at 60 mK. The quadratic dependence of dissipa- 
tion suggests the presence of charged defects or impurities 
in the samples. The scattering rate of the charged impu- 
rities is adjusted by the magnetic field and can be shown 
to give quadratic dependence j2^. Below 2 Tesla, dissi- 
pation saturates with respect to the field in both samples. 
The following temperature measurements are done at the 
saturation field of 2 Tesla. The field dependence of dissi- 
pation was measured at different temperatures up to 144 
mK and showed the same quadratic behavior. 

Table 1 summarizes calculation results found in liter- 
ature for various clamping conditions and thermoclastic 
damping. Clamping loss is due to the strain energy ra- 
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TABLE I: Quality factor from the different damping mecha- 
nisms. Here L = 6 /im, w = 300 nm and t = 200 nm are the 
length, width and tickness of the beam and h — 675 «m is 
the thickness of the substrate. 



Quality factor 



Clamping ^5J 


Lh 2 /0.95wt 2 = 2.5 x 10 8 


Clamping[25.1 h > t 


L 5 /0.31w/i 4 = 5 x 10 7 


Constriction [24J 


L/w = 20 


Thermoelastic[26j at 10 K 


10 10 



diating into the support structure of the suspended res- 
onator. Cross and Lifshitz [24| have estimated the elastic 
wave transmission across an abrupt junction in a one di- 
mensional scalar model. Their assumption of irreversible 
energy loss into the supporting pads results in a overes- 
timate of dissipation Q^ 1 = w/L. Photiadis and Judge 
[25| extended the analysis of the doubly-clamped beam 
and calculated the clamping loss. Using our beam param- 
eters, the resulting energy loss due to clamping is found 
to be on the order of 10 -8 . Thermolastic damping is due 
to the interaction of the normal modes of vibration of a 
resonator with thermally excited elastic modes. Lifshitz 
and Roukes [2(| have estimated thermoelastic damping 
for silicon in a regime where thermal phonons are diffu- 
sive. The value of Q^ 1 they report is less than 10~ 10 at 
temperatures below 10 Kelvin. The observed dissipation 
in our beams is orders of magnitude greater than both 
clamping and thermoelastic loss. Additional friction and 
self-heating arises from the metal electrode on the top 
surface of our beams, as shown in the experiments of Rcf. 
27. The latter mechanism is a nonlinear effect that be- 
comes inapplicable in our linear response regime, where 
the dissipation is independent of the driving current. We 
verify this fact by monitoring the response peak width at 
different driving voltages in the linear regime (not shown 
here) and observing no variation, within the measure- 
ment error. Furthermore, the temperature dependence 
of the data below 1 Kelvin is not expected to result from 
either of these mechanisms. We rule out thermoelastic 
and clamping loss as well as metal electrode friction as 
mechanisms responsible for the observed dissipation in 
our oscillators. 

Dissipation behavior very similar to what we observe 
(Fig. 4) has been seen for the acoustical properties and 
the dissipation in amorphous materials and disordered 
crystals (H El El H It can be described by the ex- 
istence of two- level systems (TLS), which represent de- 
fects in the crystal. The defects arise from the intrin- 
sic impurity atoms in the bond structure, the broken or 
dangling bonds on the surface due to the abrupt ter- 
mination of the crystal structure, and contamination by 
other atoms like oxygen and water molecules. Accord- 
ing to the glass model, these TLS move in asymmetric 
double-well potentials with asymmetry A and charac- 
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FIG. 3: (Color online) a) Temperature dependence of dissi- 
pation for the 12 MHz beam on a log-log scale. The T~~ 2 and 
exp(—2e/kBT)/kBT curves are the predicted acoustic dissi- 
pation forms due to TLS in crystals at low and high temper- 
atures respectively. The best fit to our dissipation data is the 
T ' 36 power law. b) Temperature dependence of the shift in 
the resonance frequency on a semi-log scale with a guide to 
the eye (dashed line). According to the acoustic dissipation 
in crystals, the temperature dependence of the frequency shift 
at higher temperatures goes as T -1 . 



TABLE II: Temperature dependence of the dissipation and 
shift in the resonance frequency arising from the glass and 
single crystal models. Here C = noD 2 / pv 2 , E is the energy 
of TLS, a = -K 3 CM 2 k 3 /24v 5 ph 4 Lo, b = 2CuD/(w% + D 2 ) and 
d — ui/k'gD. The number of TLS and the density of states 
for the TLS (J~ 1 m~' i ) for both models are given. 
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teristic tunnel splitting energies Ao on the order of the 
thermal energy, or 10~ 4 eV at IK. The TLS couple to 
the phonons of the lattice resonantly at low tempera- 
tures (typically below IK), while at higher temperatures 
the phonons modulate the asymmetry energies which re- 
sults in relaxation absorption by the TLS [ill E2I Ei^ . 
In glasses the distribution of the asymmetry energies is 
very broad and the density of states P for the TLS's is 
taken to be constant in terms of t he p arameters A and 
A , P(A, A )dAdA = j^dAdA |2g, whereas for sin- 
gle crystals Phillips has suggested that similar thermal 
properties can be obtained with a well defined tunnel 
splitting energy and a Gaussian asymmetry distribution 
of width Ai, P(A)dA = Aexp(- A 2 /2A\)dA [H. Both 
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FIG. 4: (Color online) a) The sharp response of the 12 
MHz beam as temperature is increased displays significant 
frequency shift and dissipation. The shift in the resonance fre- 
quency is illustrated by means of a guide to the eye (dashed 
line), drawn on the raw data, b) Temperature dependence 
of the shift in the resonance frequency for both beams on a 
semi-log scale. The plot for the 14.586 MHz beam has been 
shifted up for clarity. The data exhibits a peak in the vicinity 
of 0.7 K in both samples. Below the temperature 0.7 K the 
dependence is closely logarithmic InT. However, it is difficult 
to determine the behavior above 0.7 K due to the sparseness 
of the data, c) Below 2 Kelvin dissipation shows weak power 
law dependence T 0,36 . Both samples exhibit saturation of 
dissipation below 100 mK. 



approaches agree on the low temperature behavior of fre- 
quency shift and dissipation, and the dependences are 
outlined in Table 2. 

In single crystals, dissipation due to TLS is expected 
to have a T~ 2 form at high temperatures E < ksT, and 
decreases as T~ 1 e~ 2E / kBT at low temperatures E > ksT 
where E is the energy of TLS. The sound velocity is 
predicted to follow the T -1 dependence at high tem- 
peratures, saturating at lower temperatures |23l l28t l2^| . 
Fig. 3 shows poor agreement of the crystalline impurity 
dissipation theory with our data. 

In Fig. 4 we observe logarithmic temperature depen- 
dence in the frequency shift data corresponding to reso- 
nant absorption at temperatures below To = 0.7 K. The 
temperature dependence of dissipation in our data J 10 - 36 
is weaker than the expected T or T 3 dependence due to 
TLS relaxation via electron or phonon channels respec- 
tively [l5l l3fl | , or the T dependence in the Phillips model 
for a crystal with impurities. The contribution of reso- 



T<T , 



(1) 



fo To 1 

where Tq = 60 mK is a reference temperature [3lJ. C 
is defined by C = where p and v are density and 

sound velocity, Uq is the density of states for the TLS and 
D is the dcformaion potential. We have calculated the 
number of two-level systems in our silicon beams using 
D = 0.4 eV [13, p = 2330 kg/m 3 and v = 5000 m/s. We 
find that the measured dissipation is due to the coupling 
of the resonant phonon mode to as few as 50 atoms in 
the silicon beams. 

Because of its small size in the micron/nano scale and 
the associated large surface-to-volume ratio, the single- 
crystal structure appears to behave like a disordered or 
amorphous system with a wide distribution of asymme- 
try energies due to surface roughness. Although both the 
glass model as well as the Phillips model for disordcrd 
crystals account for most of the features, the low tem- 
perature saturation of dissipation (with no corresponding 
saturation of the frequency shift) cannot be explained in 
this framework. This excess dissipation at low tempera- 
tures may be understood in a more sophisticated analy- 
sis of the TLS such as the Caldeira-Leggett model in the 
sub-ohmic regime H3 ■ 

In conclusion, we have measured the dissipation and 
frequency shift in the resonance frequency of MHz-range 
nanomechanical beams at millikelvin temperatures. The 
temperature dependence exhibits partial agreement with 
both glass and disorderd crystal models of tunneling two- 
level systems (TLS), with quantum dissipation arising 
from coupling to as few as 50 TLS. A proper understand- 
ing of the quantum dissipation mechanisms in nanome- 
chanical structures at millikelvin temperatures will be 
cssencial to the realization of the quantum harmonic os- 
cillator |3 [35)]. Our results help explain the appar- 
ent upper limit on the quality factor of nano-electro- 
mechanical structures (NEMS) with decreasing system 
size. Their smaller sizes make NEMS progressively sus- 
ceptible to (quantum) dissipation by a handful of two- 
level systems arising from dangling or uncompensated 
bonds, which, even in the cleanest structures, are as un- 
avoidable as the required termination of the surface. This 
work is supported by the National Science Foundation 
(grant number DMR-0346707) under the NSF-EC Coop- 
erative Activity in Materials Research. 
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